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Abstract 


With multi-stream coannular nozzles, sometimes tones occur that may cause the nozzle to fail noise 
regulation standards. A two-stream nozzle was studied experimentally and numerically in an attempt to 
identify the sources of such tones and explore remedies. For the given nozzle configuration, sharp tones 
occurred in a range of low jet Mach numbers. The tones apparently occurred due to a coupling between 
vortex shedding from the struts, which held the nozzles and the center-body together, with various duct 
acoustic modes. A leading edge treatment of the struts is shown to eliminate the tones via disruption of 
the vortex shedding. 


Introduction 


There have been numerous previous studies addressing tones occurring with jets and shear layers and 
their interactions with structures (e.g., Refs. 1 to 5). The present investigation addresses tones sometimes 
observed with multi-stream nozzle flows. In fact, the study was prompted by persistent tones encountered 
in large-scale multi-stream nozzle tests at NASA Glenn Research Center (GRC) (Ref. 6). Such tones may 
make the nozzle fail noise regulation standards and this led to a model-scale experimental investigation. A 
0.46-scale model of a two-stream coannular nozzle, that exhibited such tones in the large-scale tests, was 
fabricated for the investigation. The effort was to understand the origin of such tones and explore possible 
remedies. This paper reports results from the model-scale experiment. 

The flow geometries of complex multi-stream nozzles involve multiple coannular ducts as well as a 
center-body. Usually, three or four struts, spaced equally in the circumferential direction, hold these 
coannular elements together. As the flow velocity is increased from zero, distinct tones may occur. With 
increasing speed a tone may disappear and new tones with different frequencies may appear. As found in 
this investigation and as explained with the results in the following, one mechanism for the origin of such 
tones is vortex shedding from the struts. The shedding apparently couples with acoustic resonance modes 
of the nozzle configuration. The result is an acoustic tone that is spectrally more concentrated than that 
from vortex shedding alone. A remedy is found. A leading edge treatment of the struts stops organized 
vortex shedding via a disruption of the two-dimensionality of the flow. This may yield a complete 
elimination of the tones. Key results from the investigation are summarized in this paper. 


Experimental Facility 


The experiments are conducted in an open jet facility at GRC. A picture of the nozzle mounted on the 
end of the facility is shown in Figure 1(a). Unless specifically stated all dimensions are given in inches. A 
30 in. diameter plenum chamber with flow conditioning units is the main component of the facility. 
Compressed air supplied from one end of the plenum discharges through the nozzle, mounted on the other 
end, into the ambient of the test chamber. The two-stream nozzle assembly is shown by the drawing in 
Figure 1(b), with annotations identifying the various components. It consists of a cylindrical inlet section 
that is attached to the plenum chamber. This section is of 4.66 in. diameter and 3.5 in. in length. The 
diameter and slope at the junction with the plenum are matched for smooth flow lines. An interested 
reader may find further description in Reference 7. 
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(d) 
Figure 1.—Experimental facility. (a) Jet rig, (b) drawing of nozzle assembly, (c) accessories, (d) inner flow 
strut with leading edge “cap”. 


With reference to Figure 1(b), the rest of the two-stream nozzle assembly is bolted to the end of the 
cylindrical inlet section. The outer nozzle has an exit diameter of 3.46 in. Compressed air from the single 
supply splits and passes through the inner and outer passages. The inner flow passes through a removable 
nozzle having exit diameter of 3.03 in. Three otherwise identical inner nozzles are used whose lip 
thickness (tj) is varied as 0.015, 0.030, and 0.050 in. Items | and 2 in Figure 1(c) represent the inner 
nozzle with the thinnest and the thickest lip thicknesses, respectively. The “removable inlet” in 
Figure 1(b) is a straight cylindrical section with rounded leading edge. Most of the data pertain to an inlet 
with length L = 2 in. Two other inlets of length L = 0.75 and 4.75 in. (items 3 and 4 in Figure 1(c)) are 
also used in order to vary the (longitudinal) duct acoustic modes and study their impact on the audible 
tones. The ratio of the flow velocities at the exit of the inner and outer channel is measured to be 
approximately unity for all hardware at all operating pressures. Items 5 and 6 in Figure 1(c) are two 
annular discs that can be placed just upstream of the struts to block either the inner flow or the outer flow. 

The inner nozzle and the center-body are held together with respect to the outer shell of the assembly 
by four struts. The struts are placed in a “crossed” shape. Thus, both the outer and the inner flows 
encounter four struts spaced equally in the circumferential direction (Figure 1(b)). A sketch of a strut 
spanning the inner channel is shown in Figure 1(d). The span (s) is 0.8 in. and the chord (c) is 0.625 in. 
while the maximum thickness (f) is 0.125 in. Also shown in this sketch is a leading edge “treatment” in 
the form of a cap attached to the strut at mid-span. The cap geometry is described further with the results 
in the following. Results from numerical analysis of acoustic waves within the nozzle, and associated 
procedure, are discussed in a separate subsection following the experimental results. 
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Experimental Results 


Sound pressure level (SPL) spectra measured at a location perpendicular to the nozzle axis, and about 
50 in. away from the exit, are shown in Figure 2. The spectral traces (0 to 100 kHz range with 1000-line 
analysis) are staggered by 10 dB and the jet Mach number (4) is indicated for each trace. A range of Mj 
from low subsonic to just over unity is covered. These data are for an inner nozzle (Figure 1(b)) with lip 
thickness of 0.030 in. The spectra display a broadband peak; e.g., at about 25 kHz at M; = 0.256 and 
90 kHz at Mj = 1.116. Its frequency continuously increases with increasing Mj. These peaks are due to 
vortex shedding from the lip of the inner nozzle, as evident from corresponding non-dimensional value of 
the frequency. The frequency (f;) of the peak is plotted in Figure 3 as Strouhal number based on the lip 
thickness and jet velocity. Also shown in this figure are corresponding data for two other inner nozzles with 
lip thicknesses of 0.015 and 0.050 in. The Strouhal number turns out to be about 0.20, a value that 
corresponds to the non-dimensional frequency of von Karman vortex shedding. The rest of the data in this 
paper pertain to the inner nozzle with 0.030 in. lip thickness. 

In addition to the broadband peak due to vortex shedding from the lip, sharp spikes are noted in the 
spectra at lower Mach numbers (Figure 2). These are accompanied by audible tones. The trends of those 
tones are shown more clearly in Figure 4 for the same nozzle configuration with spectral traces measured 
with a narrower frequency range (0 to 20 kHz with 1000-line analysis). As the jet Mach number is increased 
the frequencies of the tone vary in stages. That is, the frequency remains constant over a range of M; then 
jumps to a new value and the pattern continues. This is further demonstrated by the tone frequency data in 
Figure 5. There are four major stages as marked in the figure. With the calculated velocity within the inner 
channel near the struts (~0.388 M)) and the strut thickness, t= 0.125 in. (Figure 1(d)), the dimensionless 
frequency (Strouhal number) in each stage turns out to be about 0.3. This is not too far from the value for 
Karman vortex shedding. (The Strouhal number is often found over a range around the expected value of 
about 0.2. This is due to the fact that the effective length scale is different from what is used as well as 
boundary layer effects. In the present case, a thickness smaller than ¢ corresponding to a location closer to 
the trailing edge is likely to be the length scale, accounting for the higher value of 0.3 that was obtained by 
using #). The reader is reminded that, unlike the broadband peaks in Figure 2 due to shedding from the inner 
nozzle lip, here we are considering shedding from the trailing edge of the struts. The sharp peaks at lower 
M,, corresponding to the audible tones, are the focus of the present study. 
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Figure 2.—Sound pressure level spectra at different jet Mach number 
(Mj) for 0.030 in. lip inner nozzle. 
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Figure 3.—Strouhal number corresponding to the broadband spectral peaks 
(Figure 2) versus M,, for three different lip thicknesses of the inner nozzle 
(Figure 1(b)). 
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Figure 4.—Sound pressure level spectra at different jet Mach number for the 
0.030 in. lip inner nozzle, in a low Mj -range with expanded abscissa. 
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Figure 6.—Tone frequency versus Mj with inlets of two different 
lengths and with a case where the outer flow is blocked. 
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In an effort to understand the nature of the acoustic modes, spectral data were acquired with varying 
length of the inlet section (Figure 1(b)). Corresponding frequency data for the tones are plotted in 
Figure 6. Whereas the data in Figure 5 corresponded to an inlet length L = 2 in., data for L = 0.75 in. 
(short) and 4.75 in. (long) inlets are shown in Figure 6. Also included are data with the outer flow blocked 
using the annular disc 6 of Figure 1(c). It is clear that for any of these parametric variations the tone 
frequencies and the staging behavior basically remain unchanged. When the inner flow was blocked by 
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disc 5, Figure 1(c), the tones disappeared. Thus, the source of the tone must be in the inner flow since the 
tones occur regardless of the presence of the outer flow. Furthermore, somewhat surprisingly, it becomes 
clear that the length of the inner duct has no impact on the frequencies or the stages of the tones. The duct 
modes are addressed further with the numerical results in the following. All subsequent data in the paper 
pertain to the inlet section with length L = 2 in. 

From the discussion of the data in Figure 4 to Figure 6 it is apparent that vortex shedding from the 
struts is the likely origin of the observed tones. However, the fact that the tone frequency remains 
constant within each stage (Figure 5 and Figure 6) indicates that there is a coupling with acoustic resonant 
frequencies characteristic of the nozzle configuration. When the shedding frequency is close, a particular 
acoustic resonance mode is excited. The excitation of that mode occurs over a range of the jet Mach 
number. When the velocity is varied sufficiently the shedding frequency becomes incompatible and the 
coupling ceases. However, the frequency can now become compatible with another acoustic mode 
leading to a tone of different frequency. This explains the staging behavior. 

It should be mentioned here that various other parametric variations were examined in an effort to 
understand the origin of the sharp tones. The parameters varied included the lip thickness of the inner 
nozzle, lip-to-lip distance between the inner and the outer nozzle, a flared and a constricted inlet instead 
of the straight removable inlet in Figure 1(b), as well as the overall length of the inner duct by varying the 
length of the removable inlet as already discussed. (When designing the experiment, the origin of the 
tones was suspected to be with the flow near the nozzle exit and its interaction with the duct modes. The 
struts were not even considered as culprits and hence no provision was made for varying the strut 
geometry.) For these explorations four conditions were chosen from the four stages of Figure 5. 
Corresponding spectral traces at M; = 0.168, 0.260, 0.345, and 0.427 are shown in Figure 7, which display 
spikes at frequencies 4.5, 7.75, 9.5, and 12.38 kHz, respectively. For clarity the spikes are marked by the 
circular symbols. Variation in the lip-to-lip distance was achieved by unscrewing the inner nozzle (that, of 
course, left a cavity on the inner wall of the inner duct). None of these variations affected the two spectral 
peaks at the lowest frequencies (4.5 and 7.75 kHz, Figure 7). With lip-to-lip distance variation the spectral 
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Figure 7.—Sound pressure level spectra at four jet Mach numbers for detailed 
parametric study. 
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peaks at higher Mj; were altered with appearance of peaks at other frequencies, probably due to additional 
tones originating from the inner wall cavity. In general, the tones remained unchanged for variation of the 
aforementioned parameters and eventually it became apparent that the tones were tied to the struts in the 
inner duct. 

A key result of the study was obtained when the leading edge cap was applied on the inner struts. The 
tones were completely eliminated as demonstrated in Figure 8. The caps were cut from a plastic “edge 
protector” for car doors. With reference to Figure 1(d), the dimensions were approximately, / = 0.31 in., 

h = 0.26 in., and w = 0.64 in. These were placed centrally on each of the four struts in the inner flow. 
Experiments showed that placing the caps on the struts in the outer annular flow had no effect on the 
tones. The tones were suppressed only when the caps were placed on the inner struts. Furthermore, it was 
found that caps with much smaller span w = 0.3 in. or even 0.1 in. were just as effective in eliminating the 
tones. Apparently, the caps disrupted the two-dimensionality of the flow over the struts and thereby the 
organized vortex shedding was stopped which in turn led to the elimination of the tones. 

A further interesting result was obtained when caps with full span, equal to the span of the struts 
(w = 0.8 in.), were tried. With the full-span caps on all four inner struts the tones reappeared albeit at 
higher Mach numbers. This is shown in Figure 9. A tone at 4.13 kHz appeared at the highest M;. At even 
higher values of M; two stages of the tone could be detected, as can be seen from the spectra traces of 
Figure 10. With the no-cap condition tones were not heard at Mj; values higher than about 0.43. With the 
full-caps the occurrence of the tones shifted to a higher M; range and could be heard up to about M; = 
0.85. With the height 4 = 0.26 in. of the cap (Figure 1(d)) and the local velocity in the duct (~0.388 M4), 
the frequencies of the two stages are found to correspond to a Strouhal number of about 0.22. Therefore, 
it is apparent that with the full-span caps two-dimensionality of the flow is restored that leads to 
organized vortex shedding which in turn excites the tones. However, since the shedding frequencies are 
lower due to the larger thickness, for a given acoustic mode the frequency matching occurs at a higher 
Mach number. This is why the occurrence of the tones has shifted to a range of higher M,. 
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Figure 8.—Sound pressure level spectra at same four Mj of Figure 7 but with 
semi-span caps on the struts in the inner flow. 
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Figure 9.—Sound pressure level spectra at same four Mj of Figure 7 with 
full-span caps on the struts in the inner flow. 
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Figure 10.—Sound pressure level spectra with full-span caps on the struts 
in the inner flow, at higher values of Mj. 
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Figure 11.—Tone frequency versus Mj with full-span caps Figure 12.—Schlieren pictures of the flow; 
on inner struts. full-span caps on four inner struts. (a) Mj = 


0.45, f= 4.13 kHz, (b) Mj = 0.67, f= 6 kHz. 


For the full-caps the frequency data are plotted in Figure 11. Also shown in this figure are the 
amplitudes of the spectral peaks (in arbitrary scale) to provide an idea about the relative variation of the tone 
intensity within each stage. Clearly, the amplitude becomes the largest in the middle of each stage. Also, 
these amplitudes at higher M; with the full-caps turned out to be much larger than those encountered with 
the no-cap case. 

The higher Mach number conditions with the full caps allowed good visualization of the flow using 
the schlieren technique. Typical visualization pictures for the two stages of Figure 11 are shown in 
Figure 12. The large amplitude tones clearly excite the outer shear layer as evident from the organized 
vortices marked by the downward pointing arrows. Vortex shedding from the lip of the inner nozzle 
(leading to the broadband peaks, Figure 2) can be barely discerned as marked by the upward pointing 
arrows. Unfortunately, neither these pictures nor the experimental results presented so far provide a clear 
idea about the nature of the acoustic modes within the nozzle that are excited by vortex shedding from the 
struts. Some insight is gained by the numerical simulation results described next. 


Numerical Simulation 


The intent of this analysis is to identify duct modes that couple with Strouhal shedding from the 
struts. A commercial software package, COMSOL Multiphysics® (Ref. 8), is used to perform the analysis 
for the 2-stream nozzle configuration. A Computer Aided Design (CAD) geometry (Figure 1(b)) is used 
to create a 3D COMSOL model of the nozzle. The computed domain (Figure 13) starts upstream of the 
internal flow split between the core- and bypass-flow paths, and extends downstream to just past the tip of 
the center plug. The simulation solves the acoustic pressure physics using a finite element formulation on 
an unstructured mesh. There is no flow in the domain. The upstream boundary has a plane wave radiation 
condition. The downstream radial and axial surfaces use a spherical wave radiation boundary condition 
with center on the jet axis and axially at the inner nozzle lip. 
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Figure 13.—Computational domain. Monitored plane 
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Figure 14.—Maximum of response function within 
computed domain versus frequency of perturbation. 


identified is the driven strut. The coordinate labels 
in this figure are in meters. 


(a) Maximum pressure, (b) maximum r.m.s. velocity. 


Two small panels on either side of one of the struts in the inner flow are used to provide periodic 
perturbation. The panels are placed near the trailing edge. The “driven strut” at the top, with the panels 
marked blue, is identified in Figure 13. The two sides are driven 180 degrees out of phase in order to 
simulate the asymmetric vortex shedding. The perturbation pressure amplitude is set arbitrarily at 
+1 Pascal. 

With the perturbation at a given frequency the simulation is performed and the acoustic field in the 
domain is then interrogated. For each frequency, two metrics for resonance strength are collected: the 
maximum of the pressure fluctuation amplitude in the domain and the maximum of the root mean square 
acoustic velocity in the domain. A sweep in frequency is then conducted in steps, e.g., 25 Hz, 100 Hz, etc., 
with larger steps being used at higher perturbation frequencies. This way a spectrum of the response 
function is constructed that reveals the resonance conditions. Then, for specific resonant frequencies the 
mode shape, i.e., the nature of the pressure distribution in the domain is studied. As a first step, the mode 
shape over a cross-sectional plane just downstream of the strut (marked red and identified in Figure 13) 
is examined. 

The response functions (maximum pressure amplitude within the domain and corresponding maximum 
r.m.s. velocity) versus the perturbation frequency are shown in Figure 14. It immediately becomes clear that 
the simulation approximately captures the tone frequencies as noted in the experiment. The three peaks seen 
in Figure 14(b) roughly correspond to three of the peaks seen in Figure 7 at 4.46, 7.76, and 12.37 kHz. The 
peak at 9.76 kHz (Figure 7) is missed, however, there is a hint of additional energy around that frequency in 
Figure 14(b). Furthermore, the identified resonant frequencies are encouraging in that they exhibit 
amplification at the non-driven struts, indicating a synchronized Strouhal shedding. As discussed further in 
the following, the frequency agreement with the experiment is not exact, however, quite satisfactory 
considering the facts that the simulation is for no-flow and there are experimental uncertainties 
(approximately 100 Hz) due to variations in ambient conditions from day to day. 

The pressure distributions on the monitored plane (Figure 13) are shown in Figure 15 for the three 
dominant resonant frequencies of Figure 14(b). These frequencies are 4.525, 8.05, and 12.52 kHz in 
Figure 15(a) to (c), respectively. Each distribution represents the amplitude in phase with the driven 
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surfaces. Clearly, “standing wave” patterns have set up among the four struts. A region of either negative 
or positive pressure exists between a given pair of struts in Figure 15(a) and the pattern is anti-symmetric 
across a diameter containing the opposite struts. It appears that the azimuthal spacing of the struts plays a 
role in dictating the resonance condition. The pattern in Figure 15(a) corresponds to the tone at the lowest 
frequency stage (Figure 5) and thus may be considered as the “fundamental” in the phenomenon. 

The data for resonance at 8.05 kHz (Figure 15(b)) exhibit pairs of positive and negative pressure 
regions within an intra-strut space. Thus, this may be considered as the first harmonic in the phenomenon. 
The distribution in Figure 15(c) is also similar but there is an additional radial distribution. This 
represents a complex higher harmonic. Resonances at many such harmonics are possible. In fact, an 
eigenmode analysis using the COMSOL solver identified great many eigenfrequencies. For example, 

98 modes were found for the given geometry in the 3500 to 5500 Hz band. Only a selected few of these 
would be compatible with, and thus excitable by, the asymmetric perturbation from the TE of one of the 
inner struts. Note that due to prohibitive computational time at high frequencies only one-half of the 
cross-section is computed for the case in Figure 15(c). Furthermore, the amplified scale in the latter case 
(+1.5 Pa, compared to +3.5 Pa in Figure 15(b)) made the pressure distribution within the outer annular 
passage also visible. 


(c) 
Figure 15.—Cross-sectional distribution of pressure 
at the monitored plane (Figure 13). Frequency and 


amplitude range are: (a) 4.525 kHz, +2.03 Pa, 
(b) 8.05 kHz, +3.5 Pa, and (c) 12.52 kHz, +1.5 Pa. 
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A corresponding pressure distribution for the fundamental (4.525 kHz) over the entire domain is shown 
by six iso-surfaces in Figure 16(a), for pressures of +1, +2, and +3 Pa. Another perspective of the same 
domain is shown in Figure 16(b), skipping the +2 Pa iso-surfaces to avoid clutter. These distributions 
indicate a complex mode shape. It is not a simple longitudinal, helical or radial duct mode. From the +3 Pa 
iso-surfaces, marked in Figure 16(b), it is apparent that regions of high pressure have formed, against the 
outer wall of the inner channel, in between pairs of struts. (The walls of the inner and outer channels are 
made transparent in the figures in order to make the iso-surfaces visible.) The +1 Pa iso-surfaces (with 
lighter red and blue colors, Figure 16(b)) touch the “driven” strut since the driving pressure is of that 
amplitude. Furthermore, the latter iso-surface touches all the struts demonstrating the synchronized motion 
in the domain. There must be acoustic feedback among the struts in order to achieve the synchronization of 
the vortex shedding. Note once again that there is high pressure buildup in the intra-strut spaces. These are 
locations of pressure anti-nodes. Therefore, the struts themselves should be locations of pressure nodes and 
hence anti-nodes for induced velocity. This is indeed found to be the case as demonstrated in Figure 16(c). 
High induced velocity near the TE of the struts would be consistent with synchronized vortex shedding in 
the case with flow. 


Figure 16.—Iso-surfaces of pressure and velocity for 
f = 4.525 kHz over entire domain. (a) Side view with 
iso-surfaces of pressure at +1, +2 and +3 Pa. (b) Oblique 
view with iso-surfaces of +1 and +3 Pa. (c) Oblique view 
with an iso-surface of r.m.s. velocity at 0.006 m/s. 
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TABLE 1.—COMPARISON OF PREDICTED FREQUENCY WITH TONES OBSERVED IN THE EXPERIMENT 


Mj, Jet Mach Moocal, inner Correction, Freq, f (Hz) Freq, f (Hz) Freq, f (Hz) Difference 
number channel M at Sqrt(1- Miocai) predicted corrected expt from expt. 
strut location 
0.168 0.065 0.998 4575 4565 4460 2.4% 
0.260 0.101 0.995 8095 8054 7760 3.8% 
0.427 0.166 0.987 12690 12522 12375 1.2% 


The high pressure buildup against the inner path outer wall in the intra-strut spaces also indicate 
reinforcement by reflected waves from the wall, and hence a necessary role of the duct surrounding the 
struts to establish the standing waves. In other words, such resonance would not be expected if the walls 
were removed and the strut configuration were left in free space. Except for imposition of a cut-on 
frequency for the waves to emanate through the nozzle exit, any additional role of the ducts (such as, 
possible feedback between the struts and the nozzle exit) in establishing the resonances remains unclear at 
this time. In any case, the standing wave pattern around the struts under the artificial perturbation appears to 
be the source from which waves propagate in both upstream and downstream directions. In practice, with 
the flow on, the organized vortex shedding substitutes for the artificial perturbation to generate the tones. 

Finally, it should be noted that the resonance frequencies found from the simulation were somewhat 
higher than those seen in the experiment. It is possible this is due to the fact that the simulation involved 
no-flow. With flow, the feedback from one strut to the next would require longer wavelength, directed 
upstream at an angle, whose component on the plane of the struts would excite and synchronize with 
shedding from the adjacent struts. Thus, a simple geometric correction to account for the longer 
wavelength due to local flow velocity would yield the correct frequencies. Table 1 shows such a 
correction applied to the predicted frequencies from the no-flow simulation and the resulting agreement 
with the experimental data. 


Conclusions 


A two-stream nozzle was studied experimentally and numerically in an attempt to identify the sources 
of resonant tones. With the given nozzle configuration, sharp tones occurred in a range of low jet Mach 
numbers. These appeared due to a coupling between vortex shedding from the struts, which held the 
nozzles and the center-body together, and acoustic resonance modes characteristic of the nozzle 
configuration. With the present nozzle geometry the source of the tone is determined to be the four struts 
within the inner channel. A standing wave pattern occurred around those struts whose shape varied and 
became more complex with increasing frequency. The “fundamental” involved alternate positive and 
negative pressure cells in the intra-strut spaces, the pattern being anti-symmetric about a diametral plane 
containing the opposite struts. At higher frequencies, more complex mode shapes occurred involving 
multiple pressure cells in a given intra-strut space as well as multiple cells occurring in the radial 
direction. A leading edge treatment of the struts in the inner channel is shown to eliminate the tones. This 
occurred due to a disruption of the organized vortex shedding. Apart from the nozzle flows considered in 
this study, conceivably a similar leading edge treatment could be applicable in many other situations, e.g., 
in wind tunnel tests, where tones induced by similar vortex shedding might be a problem. The remedy is 
simple and can be included in the original design of the hardware or retrofitted with existing hardware. 
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